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.20.  ABSTRACT  (ContinuM) 

^bewn  eneigy.  We  preeent  e  genetal  telf-coniieteAt  non-Unetr  theory  of  the  free  electron  Ueer 
proceei.  The  nonlinear  formulation  of  the  temporal  iteady  ctate  ftae  electron  laser  problem  reeulti 
tn  r,  set  of  coupled  differential  equatiohi  govern^  the  vetial  evolution  of  the  amplitudea  and 
wavelength  the  radiation  and  ipace  charge  fk^lds.  ThM  equationa  are  readily  aohred  numerically 
<ince  the  amplitude  and  waviden^  very  on  a  spatial  scale  whkh  is  comparable  to  a  growth  length 
of  the  output  radiation,  A  number  of  numerical/analytical  illustrations  are  presented  ranging  from 
the  optical  to  the  stibmilUmeter  wavelenath  regime.  Our  non'Unear  formulation  in  the  linear  regime 
is  compared  with  linear  theory  and  agreement  is  found  to  be  excellent..  Analytical  expressions  for 
the  saturated  efficiency  and  radiation  amplitude  are  also  shown  to  be  in  very  good  agreement  with 
our  non-linear  numerical  solutions.  Efflciency  curves  are  obtained  for  both  the  optical  and 
submiillimeter  PEL  examples  with  flxed  magnetic  pump  parameters.  We  show  that  these  intrinsic 
efficiencies  can  be  greatly  enhaitced  by  appropriately  contouring  the  magnetic  pump  periodc^^ 
the  case  of  the  optical  PEL,  the  theoretical  s^e  pass  efficiency  can  be  greater  than  20%  b^"***^ 
appropriately  decreasing  the  pump  period  and  increasing  the  pump  magnetic  field. 
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THE  NON-LINEAR  THEORY  OF  FREE  ELECTRON  LASERS  AND 
EFFICIENCY  ENHANCEMENT 


1.  INTRODUCTION 

Free  electron  lasers  (FELs)  based  on  backscattering  from  relativistic  electron  beams  have 
demonstrated  a  unique  potential  for  becoming  a  new  type  of  coherent  radiation  source.  In  prin¬ 
ciple,  these  radiatl,-'n  sources  will  be  characterized  by  output  wavelengths  ranging  from  the  mil- 
limeter  to  beyond  the  optical  regime,  frequency  tunability,  very  high  power  levels  and  high 
eti'iciencies. 


Theoretical  analysis  on  the  FEL  mechanism  has  been  carried  out  in  the  single  particle’''* 
as  well  as  the  collective  scattering  regime.^’ Also,  non-linear  processes  and  saturation 
efficiencies  have  been  considered  for  various  FEL  scattering  regimes.*  ®  '^”  ’*  ”’^^'^’ 

The  operative  mechanism  in  FELs  is  a  parametric  process  in  which  a  long  wavelength 
pump  held  interacts  with  a  beam  of  relativistic  electrons.  Under  certain  conditions  the  incident 
pump  field  will  decay  into  a  longitudinal  wave  (density  wave)  and  a  backscattered  electromag¬ 
netic  wave  which  is  double  doppler  upshifted  in  frequency.  The  longitudinal  wave  (also  re¬ 
ferred  to  as  density  wave,  beat  wave  or  ponderomotive  wave)  results  from  the  coupling  of  the 
pump  field  and  the  electromagnetic  field  through  the  v  x  B/c  force  term.  The  ponderomotive 
wave  plays  a  central  role  in  the  linear  and  non-linear  development  of  the  scattering  process.  Its 
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effects  on  the  electron  beam  is  closely  analogous  to  the  role  played  by  the  negative  energy 
(slow  space  charge)  wave  in  conventional  traveling  vmve  mechanisms. 

T  he  pump  field  may  take  the  form  of  a  static  spatially  periodic  magnetic  or  electric  field  or 
a  propagating  electromagnetic  wave.  In  this  paper  we  take  the  pump  to  be  a  static,  periodic 
right-handed,  helically  polarized,  magnetic  field.  The  frequency  of  the  scattered  radiation  is 
given  by 

6»  s=  (1  -f.  vjc)  =  Any  t  cl  I, 

where 

y,  -  (1  - 

V.  is  the  axial  beam  velocity  and  /  is  the  pump  period.  The  possibility  of  using  a  two-stage  PEL 
scattering  process,  in  order  to  reduce  the  electron  energy  required  for  very  short  output 
wavelengths,  h'':  been  suggested.'^  '* 

Roughly  speaking,  FELs  can  be  divided  into  two  categories,  depending  on  the  gain  of  the 
radiation  field.  In  the  low  gain  regime,  the  overall  spatially  integrated  gain  is  due  to  wave  in¬ 
terference  effects  and  is  much  less  than  unity.  This  is  a  single  particle  (collective  effects  are 
not  manifested  through  space  charge  fields)  scattering  regime  and  is  exemplified  by  experiments 
at  Stanford  University.^*' 

The  high  gain  FELs  are  characterized  by  stimulated  radiation  fields  which  grow  exponen¬ 
tially  in  the  linear  regime.  Experiments  with  intense  relativistic  electron  beams  performed  at 
NRL,  Columbia  University  and  Cornell  University  fall  into  this  class.^^^  For  a  detailed 
theoretical  discussion  of  the  various  FEL  mechanisms  the  reader  is  referred  to  Refs.  (11),  (17) 
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The  mein  ottiecUvet  of  thli  work  are  to  preaeot  a  aelf*conalatent  ooa*linear  formulation  of 
the  FEL  meclumiam  and  to  theoretically  anatyxa  aome  of  the  concepti  neoeaaary  to  develop 
efflcient,  high  power,  tunaMe  PEL  radiation  aouroea.  Some  of  the  aalient  features  of  this  theory 
include:  i)  completely  arbitrary  magnetic  pump  fletd  (period  and  amplit  jde  can  be  functions  of 
axla!  position)  ii)  space  charge  effecta,  ill)  arbitrary  polariiation  of  the  radiation  field,  iv)  com¬ 
pletely  relativistic  particle  dynamics  and  v)  frequency  and  spatial  harmonics  in  the  excited 
fields.  The  non-linear  formalism  developed  for  the  FEL  problem  is  also  applicable  to  a  large 
class  of  temporal  steady  state  convective  processes.  Our  formulation  of  the  problem  permits 
the  rpatial  dependence  of  the  pump  magnetic  field  to  be  arbitrary.  Hence,  efficiency  enhance¬ 
ment  schemes  which  utilize  amplitude  and  wavelength  spatial  variations  of  the  pump  field  can 
be  analyzed.  The  spatial  variation  of  the  scattered  radiation  amplitude  and  wavelength  occurs 
on  a  scale-length  which  is  large  compared  to  the  wavelength  of  the  pump  field.  This  permits 
numerical  solutions  for  cases  where  the  electron  beam  energy  is  extremely  high.  That  is,  in 
this  approach,  there  is  no  large  separation  of  spatial  scale  lengths,  despite  the  large  sp'\tial  scale 
difference  between  the  wavelength  of  the  scattered  field  and  the  pump  field,  so  arbitrarily  high 
values  of  the  relativistic  gamma  factor,  y,  associated  with  the  beam,  can  be  considered.  Fur¬ 
thermore,  the  formulation  is  carried  out  in  the  laboratory  frame  under  temporal  steady  state 
conditions. 

The  analytical  formulation  of  the  general  non-linear  steady  state  FEL  problem  consists 
essentially  of  three  parts.  In  section  II,  the  wave  equations  are  used  to  derive  expressions  for 
the  slow  spatial  evolution  of  the  amplitudes  and  phases  of  the  scatterec  fields  in  terms  of  the 
driving  currents.  Then,  in  section  III,  the  driving  currents  are  expressed  as  functions  of  the 
dynamics  of  the  particle  ensemble  (electron  distribution  function).  The  particle  orbit  equations 
are  written  selfconsistently  in  terms  of  the  scattered  fields  in  section  IV.  The  orbit  equations 
describing  the  motion  transverse  and  parallel  to  the  electron  stream  are  completely  decoupled. 
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The  Hnoer  spatial  growth  rates,  efficiencies  and  saturated  field  amplitudes  are  derived  in  section 
V  Tor  various  scattering  regimes.  Finally  a  iiumber  of  analyticui/nuraerical  illustrations  in  the 
high  gain  scattering  regime  are  given  in  section  VI.  The  non*linesr  particle  dynamics  is  dis* 
cussed  in  some  detail.  Efficiency  curves  are  obtained,  atid  a  method  of  dramatically  increasing 
the  single  pass  efficiency,  as  suggested  in  Ref.  <17),  by  contouring  both  the  pump  period  and 
magnetic  fluid  is  analyzed.^*  For  instance,  efficiency  of  >20%  are  shown  to  be  theoretically  pos* 
sible  at  optical  wavelengths  using  this  approa9h.  The  basic  idea  is  to  grudually  slow  down  the 
phase  velocity  of  the  ponderomotive  wave  at  the  point  where  the  electrons  are  deeply  trapped 
in  the  ponderomotive  wave  potentials.  The  slowing  down  of  the  wave  is  accomplished  by  adia- 
batically  decreasing  the  pump  field  period.  Appendix  A  contains  the  formulation  of  the  FEL 
process  with  spatial  and  temporal  harmonics  in  both  the  radiation  and  space  charge  fields.  Also 
included  in  this  formulation  is  the  ability  of  the  radiation  field  to  undergo  a  change  in  polariza¬ 
tion  from  a  circularly  polarized  to  an  elliptically  polarized  wave  as  the  particle  dynamics  become 
non-linear. 

The  resulting  set  of  non-linear  coupled  equations  self-consistently  relate  the  spatial 
dynamics  of  the  particles  and  fields.  These  equations  have  been  solved  analytically  in  the  linear 
approximation  and  the  linear  dispersion  relation  was  obtained.  The  full  set  of  non-linear  cou¬ 
pled  equations  are  readily  solved  numerically  for  the  spatial  growth  rate  and  saturation  level  of 
the  scattered  fields. 

II.  NON-LINEAR  EVOLUTION  OF  SCATTERED  WAVES 

The  physical  mode!  which  we  will  develop  is  that  of  a  fully  relativistic  electron  beam  in¬ 
teracting  with  a  spatially  periodic  pump  magnetic  field  as  depicted  in  Fig.  (1).  Only  spatial  vari¬ 
ations  along  the  z  axis  will  be  considered  for  the  electron  beam,  pump  field  and  scattered  radia- 
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The  variable  aiuptitude  and  period  pump  magnetic  field  can  be  expressed  in  terms  of  the 
vector  potential 

A.(r)  -  jcos  (/J  /t.(t')dr*j  -h  sinj/^'  *.(2  )*')  A  (1) 


where  the  amplitude  A^{z)  and  wavenuniber  k^{z)  are  known  and  slowly  varying  functions  of 
2.  The  potential  field  in  (1)  is  a  good  approximation  to  a  right*handed  polarized  helical  magnet¬ 
ic  field  near  the  2>axis  of  an  appropriate  coil  winding.  The  pump  magnetic  field  associated  with 
(1)  is  given  by 


B,(2) 


where 


BAz) 


\s: 


cos  1  L  k,(2')<fe'  + 


'^(2)] 


(2) 


BAt)  -  -  ((*,(2)>4c(*))^  +  (d^,(2)/d2)*)»'* 


and 


^(r) 


,  f  di4o(2)/62 
-  -  ran  '  1 — - -r- 


k„(2)  ilo(2) 


are  slowly  varying  functions  of  2.  The  period  of  the  magnetic  field  is 


1(2)  -  2n7(A<,(z)  +  BvIBz)  =  Infkoiz)- 


(3) 


The  scattered  electromagnetic  and  electrostatic  fieids  in  terms  of  the  vector  potential 
A(z,/)  and  scalar  potential  ^(z,r)  are  taker  to  be 

A(z,0  AAz)  cosjj^  /c+tz')dr’  —  o»f  +  oj  —  .<4v(r)sin|j^  k+(z')<fe'  —  wf  +  fl)|  e^,. 

^(z,»)  -  ^(z)  cosjj^  cCt')*  -  vi  +  e,j.  (4a, b) 

where  the  amplitudes  of  the  potentials.  snd  dCz)  as  well  as  the  wavenuntbers 

kA^)  snd  kU)  are  slowly  varying  ftinctioos  of  z.  'iTie  scattered  electromagnetic  field 
represented  by  Eq.  (4a)  is  a  right-handed  elliptically  polarized  field  traveling  towards  the  right. 
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X  » 

The  frequency  «  of  the  field  end  the  pheees,  f  and  ere  independent  of  r.  In  appendix  A,  the 
general  form  for  the  scattered  field!  la  used  in  the  non*Iinear  formulation  of  the  problem. 

However,  in  the  main  body  of  the  text,  the  fields  in  Eqs.  (4)  will  be  used  in  order  to  minimise 

the  notational  algebra.  ,  i 

{ 

The  evolution  of  the  scattered  potentials  is  governed  by  the  wave  equations 

A  (s.f)  -  -  ! 

j 

and 

“  ‘♦**'.<**')*  <5a.b) 

dadf 

where  J(z,r)  is  the  driving  current  density.  Substituting  the  potentials  of  Eqs.  (4)  into  (S)  we 
obtain 


_ L  ^ 

fix*  c*  Sr* 


(«Vc*  -  **  (x))  /lx(x)  cos  ♦(x.r) 

(x) -I-  (^,(x)*y* (x))  sin  ♦(x.r)  -  -  —  (x.r) , 

ox  c 

(juVe*  -  *|(x))  Af(z)  sin 

A  (i4y(x)  *y*(x))  cos  ♦(x,r)  ^  —  Jy  (z,t), 

ox  C 

sin  tlfj(z,t)  +  k(z)  ^(x)  cos  ^j(z.t)  —  —  y,(x,r),  (6a,b,c) 

02  ^ 

where 

•Kz.t)  -  k+(z')  dz'  -  oir  +  9. 

and 

♦x(x,r)  -  ir(x')  dz'  ~mt  + 

Terms  proportional  to  have  bee.*',  neglected  from  (6).  The  neglect  of  b^A/Bz^  terms  is 

not  central  to  our  formulation,  though  it  can  be  shown  to  be  an  excellent  approximation  which 
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timpliflM  tho  final  FBL  oquationa.  The  coefflcianu  of  the  linuaoidal  terma  on  the  left  hand 
aide  of  Bqa.  (6)  are  alowly  varying  ftincUona  of  t  and  independent  of  t  The  argumenta  of  the 
aintiaoMai  terma  on  the  other  hand  are  rapidly  varying  ftinetion  of  t  for  t  fixed.  The  rapidly 

Goa  ,  ^ 

time  varying  terma,  in  for  example  Bq.  (6a),  can  be  removed  by  multiplying  it  by 

^  ifi.  Performing  thia 

operation  on  Eq.  (6a)  aa  well  aa  aimilar  onea  on  B<^.  (6b,c)  we  obtain 

(wVc*  -  klit))  AAt)  -  y,  (fti)  coa  ♦(*.#)  dr. 

2*i'*(2)-J-  (yl,(x) *!'*(*))  -  ~  J,  (x.r)  ain  ^(x./)  dt, 
ox  C  *'«> 

(«Vc*  -  kl  (x))  A,  (x)  -  ~  Jy  (x,  i)  ain  ♦(x.r)  dr. 

2*y^(x)  (AyU)  klf^(x))  -  ^  //''*  y,(x.r)  cos  ♦(x.r)  dt. 

-  4  y.(x.r)  ain  ♦,(x.r)  dr. 

^  Jf(z.t)  COS  ♦,(x.r)  dr.  (7a-0 

III.  DERIVATION  OF  NON-LINEAR  DRIVING  CURRENTS 

It  is  now  necessary  to  derive  expresricns  for  the  x,  y  and  x  components  of  the  current 
densities  and  perform  the  time  integration  specified  in  Eqs.  (7).  In  general  the  non-thermal 
electron  distribution  function,  written  in  terms  of  the  electron  orbits,  is 

fit.  p.  r)  -  fiQ  v,o  *(x  -  f (ro.  r)  E  (p,  - 1»,  (ro.  r)) 

kip,  -  riyit^.  i))  9ipt~  7>,(ro,  r))dfo  (8) 

where  mq  is  the  uniform  particle  density  to  the  left  of  the  interaction  region,  i.e .  x  <  0,  v^o  is 
the  constant  axial  electron  velocity  for  z  <  0,  f  (ro,  r)  is  the  axial  position  of  the  particle  at 
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time  /  which  crossed  the  z  *  0  pline  at  time  fo  ^(/o.  t)  is  the  momentum  vector  of  the 
particle  at  time  t  which  crossed  the  z  *  0  plane  at  time  /q-  Thermal  effects  which  are  charac* 
teristic  of  actual  electron  beants  can  be  easily  included  by  appropriately  modifying  the  electron 
distribution  function  in  (8).  The  integral  over  to  in  Eq.  (8)  takes  into  a^unt  the  continuous 
flow  of  particles  into  the  interaction  region.  The  current  density  associated  with  this  electron 
distribution  is 


J(r.  /) 


/(z.  p.  t)d^p 


(9) 


where  yiy)  -  (1  As  vrill  be  seen  later  it  is  necessary  to  rewrite  Eq.  (9)  in  the 

form 


where 


J(z.  /) 


-kl  hqv.q  jiito.  r)S(f  -  T(fo,  z)) 


(10) 


r(to.  2)  ■ 'o  +  X  “TT^^ 

Jd  VjCto.  z  ) 

is  the  time  it  takes  a  particle  to  reach  the  position  z  if  it  entered  the  interaction  region,  z  0,  at 
time  to  and  v.(io,  z)  is  the  axial  velocity  of  a  particle  at  position  z  which  was  at  z  ~  0  at  time  to- 


The  quantity  d(  (tg,  i)/9t  is  the  axial  velocity  v,  of  a  particle  at  time  t  which  was  at  z  -  0 
at  time  ig.  Clearly,  for  J(z,  1)  to  be  finite,  should  not  vanish  in  the  interaction  region.  If  w, 
vanishes  and  particles  are  turned  around,  multi-streaming  develops  and  the  entire  concept  of 
exp  (-/'&»/)  being  the  only  time  dependence  is  undermined  (due  to,  for  instance,  two-stream 
instabilities).  We  assume  here  that  no  particle  is  slowed  down  to  zero  velocity  in  the  laboratory 
frame,  hence 


ytTjC/o.  f))  molOf(fo.  r)/d/|  -  '»>.(fo.  t). 


(12) 
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Substituting  (12)  into  (10),  the  general  form  for  the  driving  current  becomes 
J(z.  /)  -  -  |t  I  no  »’70  ~ 

Substituting  the  above  form  for  J(z.  t)  into  the  right-hand-side  of  Eqs.  (7),  we  obtain  the  self- 
consistent  amplitudes  and  phases  of  the  scattered  potentials  in  terms  of  driving  currents.  To 
show  how  the  right-hand-sides  of  Eqs.  (7)  can  be  reduced  by  using  Eq^  (13),  we  simplify  Eq. 
(7a)  as  an  illustration. 


Substituting  the  x-component  of  Eq.  (13)  into  (7a)  gives 

J»2n/ot 

<//  J_^  (ito  Gx  (to.  2.  t)  8  (/  -  rito,  z)),  (14a) 

where  /  v 

G,Uo>  2.  ' 

Since  the  system  of  particles  and  fields  are  in  the  temporal  steady  state,  particles  which 
cross  the  z  *  0  plane  separated  in  time  by  lnlo>  will  execute  identical  orbits  which  are  separat¬ 
ed  in  time  by  2ir/(u.  It  is.  therefore,  possible  to  define  an  initial  beam  segment,  "beamlet,"  for 
which  all  possible  steady  state  orbits  of  the  actual  beam  particles  at  e  represented  by  the  particles 
in  the  beamlet,  but  are  displaced  in  time.  The  axial  length  of  the  beamlet  is  clearly  27rv,o/c<>. 
With  these  considerations  in  mind  we  find  that  the  function  G^(/„,z,r)  and  rUg,  z)  have  certain 
periodic  properties  in  their  arguments  which  permit  the  integrals  in  (14)  to  be  greatly 
simplified.  Specifically  we  note  that 


^  /  .  ^  .  IttN  .  ,  2itN 

Gxitg,  z.t)  ■■  <7^  /„  +  — ■ — ,  z,t  + 


and 


(0 


'V‘o 


,  z)  -  Tltg  + 


InN 


fi> 


z  - 


lirN 


CD 


(15a, 0) 


where  N  0,  ±1,  ±2,  ....  The  /  integration  in  Eq.  (14a)  is  over  one  wave  period,  0  to 
2n/u.  From  Eq.  (15b)  we  see  that  over  this  range  of  t  the  argument  of  the  delta  function  will 
vanish  over  an  interval  in  tg  equal  to  2n/(i>.  Therefore,  it  is  not  necessary  to  integrate  over  tg 
from  -oo  to  +00.  Finally  from  the  property  of  GxUg,  z,t)  expressed  in  (15a)  we  find  that 
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0  dt  z,t)  60 -  Htp,  z)) 

j>  2*/m 

^  GgOo,  X,rOo,  i))  dte>  (16) 

This  can  be  seen  most  euily  in  a  diagram,  of  the  region  of  integration  in  Fig.  2,  where  the  en* 
tire  0,  to)  plane  is  broken  into  squares  of  2rr/ii>  on  a  side.  Because  of  the  symmetry  property 
expressed  in  Eqs.  (15a)  and  (ISb),  the  value  of  the  integrand  is  unchanged  along  a  diagoual; 
this  is  indicated  by  certain  squares  having  the  same  letter.  Ciearly  then,  an  integral  In  the  verti> 
cal  direction,  over  the  shaded  squares,  is  the  same  as  an  integral  in  the  horizontal  direction, 
over  the  slashed  squares.  Substitutins  (16)  together  with  (14b)  into  (14a)  results  in  a 
simplified  form  for  Eq.  (7a).  All  the  integrals  on  the  right  hand  side  of  Eqs.  (7)  can  be  re< 
duced  in  exactly  the  same  way.  Doing  this  we  find  that  Eqs.  (7)  can  be  put  into  the  form 


(«Vc^  -  k}  (z))  Ajf(z) 

,11  v,o  rOo,  z))  ,,  ,  . 

-  4  i,i ,))  ““  *.■ 

Ikl'Hz)  *y=(i)) 

,,  ,  v,o  f »*/- «j*(r„  T(t„  r))  ,  ,  ,,  . 

--4|,|,,— «X, 

(<uVc*  -  kl  (z))  AyO) 

,11  v,o  f  2*/*  VyOp.  rOp,  z))  ,  , .  .  . .  ^ 

-  -4  |e|  ^  a»  tj,(r„,  rOo,  z))'  *  "  dtp. 

IkUHz)  ^  (44,(z)  *yMz)) 

,11  »'i0  f  2*/-  %(t*.  t(/,.  z))  . 

.  -4 1«|  n,  Xj  ^  •1“  ♦,Ut(i,j))  it,. 


j4  2«/» 

cos  ll»,(z,T(ta,z))  rfr,. 


(17a-0 
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Notice  that  on  the  riRht-hand-iiide  cf  the  above  equations  the  single  integrals  over  are  trom  0 
to  3tr/w.  As  we  will  see,  these  integrals  can  be  evaluated  numerically  by  followii:!  the  orbits  of 
a  relatively  small  number  of  particles  which  enter  the  interaction  region  in  any  singie  time  in¬ 
terval  of  duration  l-alft.  tjpnn  deriving  the  general  orbit  equations  for  the  particle  ensemble, 
in  the  next  section,  we  will  assume  that  the  scattered  electromagnetic  wave  is  circularly  polar¬ 
ized,  i.e.,  ^x(z)  **  Ayiz).  This  assumption  is  clearly  not  centra!  to  our  formulation. 


IV,  PARTICLE  ORBIT  EQUATIONS 

tVe  now  express  the  particle  orbits,  which  are  needed  for  the  evaluation  of  Eqs.  (17)  in 
terms  of  the  new  independent  variables  tg  and  z.  The  forces  exerted  on  the  electrons  arise 
from  the  pump  aud  scattered  potentials  given  in  Eqs.  (1)  and  (4).  We  immediately  note  that 
the  transverse  canonical  momenta  of  the  particles  is  conserved.  Therefore,  if  both  the  pump 
and  scattered  fields  are  zero  as  z  -oo,  the  transverse  particle  momenta  are  given  by 

Pxiz,  t)  “  (s4ar(r)  +  Axii.,  t)), 

c 

and 

Pyiz,  t)  -  ^^iA^yiz)  +  Ayiz,  0).  (18a,b) 

Using  Eqs.  (18)  the  longitudinal  component  of  the  force  equation  can  be  put  into  the  form 


— -V  -f  (A,(z)  +  A(r.  r))»  -  2y(z.  t)  ^  «(z,  t)  .  (19) 

dt  lyiz,  t)  [Bz  ”  fir  ] 

where  p,iz,  t)  is  the  axial  momentum  and  the  relativistic  gamma  factor  is 

yiz,  0  -  1  +  (AJz)  +  A(z,  f))*  f  .  (20) 

Equations  (18),  (19)  and  (20)  specify  the  particle  dynamics  in  terms  of  the  pump  and  scattered 


fields.  The  transverse  and  longitudinal  particle  motion  is  formally  decoupled.  To  write  Eqs- 
(18)  and  (19)  in  terms  of  the  new  independent  variables  z  and  we  note  that 
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J.i'  ■■ 

dz'htUo,  z'). 


—  —  K  (-  t)  — 

j  fi 

where  K(z,  t)  ’»  vAl.z)-  Note  that  ■—  (which  foHows  a  particle  orbit)  ^  — -  (which  is  taketJ 

dz  az 

ai  constant  time). 


In  terms  of  zand  r  we  simply  get 


Pj(tz,  t)  -  ^  (.'^m(^)  +  >4;((z,  r)), 


/>y(z,  t)  “  (yloytz)  +  i4^(z.  t)), 


t)  _  |gP 
dz 


-^(Ao(z)  +  A(z,  t))^  -2y(z,  t)  <p(z,  t) 


»21a,b,c) 


We  have  expressed  the  particle  orbits  in  terms  of  the  entry  time  t„  and  axial  position  z.  Note 

t 

that  oiir  definition  of  the  momenta  implies  ihat,  yixUo-'^^  “  Px(-\  VyOo>  t)  —  p^(z,  t),  and 
Tj,(/o,  t)  ■=  pAz,  t)  “  y(,z.  t)  Wo  Ij(.z,  t).  At  this  pohH  we  take  the  scattered  electromagnetic 
wave  to  be  circularly  polarized  and  set  ,t,(z) '=*.4^(z)-/4  (z).  To  obtain  the  final  set  of  equa¬ 
tions  for  the  amplitude  Aiz)  and  wavenur^^ter  k.t.(z)  we  first  combine  Eqs.  (17a)  and  (17b) 
with  Eqs.  (17c)  and  (17d)  respectively.  Using  the  expressions  for  and  tjj,  given  by  Eqs. 
(21a, b)  we  arrive  at  tne  following  expressions 


(wVc^  -  kliz))  A  (z)  ^  WoV,o  f  '  i7r'(/„,  t(/o,  z)) 

IC  If 


•  Zir/w 
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?.kyH2)  kl'Hi))  -  ~  m,v,o  f  ffr'Uo,  rUo, 

|y<,(2)  sinj  J^*(fc+(2')  +  *„(2'))  <fe'  -  «»T(fo,  ?)  +  »)J  dt^.  (22a, b) 

where  we  have  used  Eqs.  (1)  and  (4a)  for  K„(z)  w'd  k{z,  t)  and  w*  -  (4ir  |cP  njnto)'^^. 

For  completeness  we  rewrite  Eqs.  (17e)  and  (17f)  tor  the  scalar  potential 

M?)  «  Z!:!  Jli®  sin  If'  k{z’)dz'  -  wT(r,,z)  +  d/.. 

hz  IT  |e|  ' 

h7)  *m  -  V  lif  “  y"’ 

The  relevant  particle  dynamics  is  contained  in  Eq.  (21b)  which  is  rewritten  jn  the  form 
dnjUo.T)  -  leP  Id  /A  .  A  _  >>^1,  .rZsL.  (24) 


where 


.  .  y(z',T(ro,r'))  nto 


(Ao(z)  +  A(z,t))2  .  ^^2(2)  +  ^2(2) 


+  2/lp(z)>t  (z)  cos  [j^*(Mz')  +  k,U'))dz'  -  WT  +  d). 


(25a,b,c) 


The  non-linear  formulation  of  the  FEL  is  fully  described  by  Eqs.  (22),  (23)  and  (24).  The 
ponderomotive  potential  piays  a  central  role  in  axially  bunching  the  electron.  From  Eq.  (24) 
we  see  that  this  potential  is  given  by 

4,^^<2.t)  -  A„(2)A(2)  cos  (//  (*4^')  +  koU'))  4  (26) 
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The  amplitude  and  phase  of  the  scattered  fields  as  well  as  the  axial  beam  momentum  all  vary 
with  a  characteristic  axial  length  which  is  much  longer  than  the  pump  wavelength  L  This  fact 
allows  for  inexpensive  numerical  simulations  to  be  performed  in  the  laboratory  frame  with  ex¬ 
tremely  high  gamma  electron  beams. 

To  see  that  the  system  quantities  vary  on  a  scale  length  long  compared  to  4  we  note  that 
the  characteristic  length,  as  estimated  from  the  arguments  of  the  sinusoidal  terms  on  the  right- 
hand-side  of  Eqs.  (22),  (23)  and  (24),  is  roughly  equal  to  L  =  (*+  +  kg  -  6»/v*)“*.  However, 
since  the  '^requency  of  the  scattered  radiation  is  w  -  cit  =  (1  +  p,)  kg,  we  find  that 

L  »  l/kg  -  //2ir. 

This  fact  permits  us  to  solve  numerically  the  PEL  equations  for  arbitrarily  high  gamma  beams. 
The  more  conventional  simulation  approaches  suffer  from  the  problem  of  large  temporal  or  spa¬ 
tial  scale  differences  even  in  the  beam  frame  of  reference. 


To  complete  our  formulation  of  FELs  we  need  an  expression  for  the  efficiency.  The 
efficiency  can  be  defined  as  the  ratio  of  the  electromagnetic  energy  flux  increase  to  the  initial 
electron  energy  flux,  that  is 

c  <E(z,0  X  B(z,r)>,-  <E(o.r)  x  B(o,f)>, 

-  T- - - } - 7T - 2 - .  (27) 

4w  v„rto(yo-l)m,c* 

where  E  —  c~'dA/df,  B  -  Cj  x  dA/dz,  <•••>,  denotes  an  average  over  the  field  period  2w/a>, 
and  Hg  and  arc  Iftc  initial  beam  axial  velocity,  density  and  total  gamma  factor. 


Using  the  vector  potential  in  Eq.  (4a)  and  taking  the  radiation  field  to  the  circularly  polar¬ 
ized,  i.e.,  Ay  ^  A,  the  efficiency  in  Eq.  (27)  takes  the  form 


V  - 


m„c 


„  (k+(z)A*(z)  -  k+(0)A*(0)) 


fjo  (y  0  i ) 


(28) 


and  is  maximum  when  the  radiation  fields  saturate. 
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V.  DERIVATION  OF  LINEAR  GROWTH  RATES,  EFFICIENCIES 
AND  SATURATION  FIELD  AMPLITUDES 

In  this  section  we  present  the  salient  features  of  the  PEL  in  the  linear  regime.  Results  for 
the  linear  growth  rate  and  expressions  for  the  saturation  efficiency  and  radiation  amplitude  are 
obtained  in  the  high  gain  case,  i.e.,  where  the  radiation  field  amplitude  has  e-folded  at  least  a 
few  times.  For  a  more  detailed  derivation  of  these  quantities  see  Refs.  (17)  and  (18)  . 

In  the  high  gain  linear  regime  the  excited  space  charge  and  vector  potentials  are  of  the 

form 

♦(2.0  - +  C.C.. 

A(2.r)  -  +  tSy)  +  C.C..  (29) 

where  ♦(O)  and  il(0)  are  the  potential  amplitudes  at  the  input  end  of  the  interaction  region, 
zxQ,  and  the  wavenumbers,  k  and  k+ are  complex  and  independent  of  2.  For  a  magnetic  pump 
of  the  form  in  Eq.(l)  with  constant  amplitude  and  period  and  cold  electron  beam,  the  disper¬ 
sion  relation  is 

D(w,  k+)  ((at-v,yk)^  - 

where  /)(«,  k)  —  o»*  -  c^k^  —  /t+  —  fc  -  kg  is  the  pump  wavenumber, 

Wft  “  (4irle|*n„/»io)'^^  is  the  beam  plasma  frequency,  v„  is  the  axial  beam  velocity, 
Voi  ~  \e\Bo/(y„mgCkg)  is  the  transverse  beam  velocity,  Bg  is  the  pump  amplitude, 
To  -  (1  -  bL  -  Boi  “  Voi/c,  /3„  -  v,g/c  and  y„  -  (1  -  Pro) Since  the  elec¬ 

tromagnetic  wave  approximately  satisfies  the  dispersion  relation  w/c  »  k+  we  can  replace 
D(tti,  k+)  and  D(u.  k)  by  -2k^{k+ -  -  <aUygV^^/c)  and  -2k  kgC^  respectively.  The 

dispersion  relation  can  now  be  put  into  the  simple  form 

6k  (6k  +  2(  kg/y„)i6k  -  A*)  -  -a^kg/2.  (31) 
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where  k  ^  u/v„ (  kjy„  tk,  ik  is  complex,  |8/tl  «  it,  f  c  U. 

A*  -  *0  -  w/(2cy^),  and  o*  -  ({/ini  *«.)*•  Equation  (31)  assumes  that  the  beam  is  relativistic 
Via  ~  c  and  a>  >  >  Two  distinct  regimes  can  be  distinguished  from  the  dispersion  re¬ 

lation  in  (31). 

a)  Weak  Magnetic  Pump  Limit 

For  a  pump  magnetic  field  strength  such  that  <<  Pmi  =  4({/y^)'^^  the  space  charge 
potential  dominates  the  ponderomotive  potential  and  collective  effects  play  an  important  role. 
In  this  regime  of  scattering  the  dispersion  relation  in  (31)  yields 

6k  -  -y  -  y  Va*  yj(  -  (A*)'  (32) 

for  the  growing  root.  Maximum  spatial  linear  growth  occurs  when  there  is  no  frequency 
mismatch,  i.e.,  Ak  -  0  and  is  given  by 

r„,„  -  -  /m  (8A:)«„  -  i  (33) 

b)  Strong  Magnetic  Pump  Limit 

In  this  regime,  defined  by  the  condition  /S^j.  >  >  /3„,„  space  charge  forces  are  dominated 
by  ponderomotive  forces.  This  is  a  single  partice  scattering  regime  and  Eq.  (31)  reduces  to 

{6k)H6k  -  A*)  -  -  o*  kj2.  (34) 

The  maximum  spatial  linear  growth  rate  according  to  (34)  occurs  for  exact  frequency  matching, 

i.e.,  Ak  ~  0  and  is  given  by 

Tm.,  -  -  /m  ^  ({/3.i)*'^k„.  (35) 

To  obtain  est:  mates  for  the  saturation  levels  in  the  high  gain  regimes  we  resort  to  heuris¬ 
tic  arguments  based  on  electron  trapping  dynamics.  It  can  be  argued  that  at  saturation,  when 
electrons  are  deeply  trapped,  the  axial  velocity  of  the  electron  beam  has  decreased  by  the 
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amount  2&v  where  Av  -  ~  is  the  difference  between  the  equilibrium  axial  beam  velocity 

and  the  initial  phase  velocity  of  the  total  lonfitudinal  wave,  i.e.,  •  ttfReik).  The  decrease 

in  the  particle  kinetic  energy  is 

ss  (36) 

so  the  ener;j  conversion  efficiency  becomes 


2yiAv/c 


-2yi(v„-w/ffe(/t))/c.  (37) 

Substituting  Re(k)  *  «/v„  +  ikjyn  +  Rei^k)  and  w  iyl^cko  into  (37),  the  expression  for 
efficiency  becomes 


n  -  +  Rei6k)/k,.  (38) 

where  Re(6k)  is  determined  from  the  solution  of  the  dispersion  relation  in  (31).  Using  Eqs. 

(31)  and  (38)  we  find  that  in  the  weak  pump  and  strong  pump  limit  the  growth  rate  maximizes 

when  Ak  >  0,  i.e.,  u  -  2y^cke,  and  the  efficiency  at  saturation  is  respectively 


and 


■n  -  €/y». 


(39) 


1,  -  +  f/y».  (40) 

Applying  the  conservation  law  for  total  energy  ffux  we  find  that  the  magnitude  cf  the  vector  po' 
tential  at  saturation  is 


where  i)  is  given  by  either  (39)  or  (40)  depending  on  whether  the  weak  or  strong  pump  limit  is 
applicable. 
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VI.  RESULTS  AND  DISCUSSION 

In  this  section  we  present  the  numerical  results  for  the  coupled  nonlinear  FEL  equations 
in  (22),  (23)  and  (24).  Illustrations  for  a  wide  range  of  parameters  ranging  from  the  submil* 
limeter  to  the  optical  radiation  regime  are  given.  The  mono-energetic  electron  beam  enters  the 
interaction  region  at  z  —  0  with  a  uniform  density.  The  magnetic  pump  field  given  in  (1)  is  as¬ 
sumed  to  be  built  up  adiabatically  from  z  <  0  to  its  initial  value  at  z  «  0.  In  all  of  our  numeri¬ 
cal  simulations  a  small  amplitude  radiation  field  is  introduced  as  a  perturbation  at  z  -  0  and  al¬ 
lowed  to  grow  spatially  and  self-consistently  accoiding  to  the  FEL  equations.  The  small  initial 
radiation  field,  typically  less  than  0.1%  of  the  saturated  field  amplitude,  allows  for  a  long  spatial 
region  of  linear  interaction  and,  hence,  for  an  accurate  comparison  with  the  linear  theory 
presented  in  section  V.  Furthenr.oiu,  apace  charge  fields  are  included  in  all  of  our  numerical 
illustrations  even  though  in  some  cases  the  ponderomotive  field  may  dominate  the  process  as  is 
the  case  in  the  strong  magnetic  pump  scattering  limit. 

We  will  first  consider  two  examples  where  the  magnetic  pump  parameters  are  fixed,  i.e., 
constant  amplitude  and  periot.  Fuither.  iore,  we  will  show  that  efficiency  can  be  increased  to  a 
few  10*s  of  percent  even  in  the  optical  regime  by  contouring  the  mag.'ietic  pump  period  .>nd 
amplitude. 

a)  Constant  Magnetic  Pump  Illustrations 

Two  examples  will  be  discu-iscvi  in  some  detail:  1)  optical  radiation  at  A  ~  0.7S  iita  from 
a  66  MeV  electron  beam  and  2)  submillimeter  radiation  at  X  338  fim  with  a  2  6  MeV  elec¬ 
tron  beam.  Table  I  lists  the  salient  parameters  for  the  magnetic  pump,  electron  beam  and  out¬ 
put  radiation  of  both  examples. 
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For  the  optical  radiation  case,  example  1,  the  matnetis  pump  amplitude  is  6.0  kO  and  the 
period  is  Axed  at  l.S  cm.  The  66  MeV  >  131),  2  kA  electron  beam  has  a  transverse  equili¬ 
brium  velocity  of  ■*  6.4  x  10'^  c  with  the  fiven  value  of  matnetic  pump  field.  The  critical 
transverse  velocity,  see  section  V,  is  •  l.S  x  10~^  c,  hence  the  scattering  process  is  in  the 
strong  pump  regime. 

Figure  (3)  shows  the  amplitude  of  the  vector  potential  of  the  scattered  radiation,  i4(z), 
and  the  spatial  growth  rate,  T  ->  0(lny<  (z))/0z,  as  a  function  of  z.  Those  plots  are  for  an  optical 
frequency  of  at  *  lyl,  ckg  *  2.S2S  x  10’^  sec~'.  Notice  that  in  Fig.  (3)  there  is  a  long  spatial 
region  where  the  growth  rate  is  fairly  constant.  This  is  the  linear  region  of  the  Interaction.  The 
value  of  the  radiation  frequency  in  this  figure  has  been  chosen  to  maximize  the  linear  growth 
rate,  i.e.,  zero  frequency  mismatch,  6k  ->  0.  The  linear  e-folding  length  associated  with  this 
output  frequency  is  38  cm. 

Figure  (4)  shows  a  comparison  between  the  spatial  growth  rates  obtained  from  the  linear 
regime  of  the  numerical  simulation  of  our  FEL  equations  (crosses  (x))  and  the  linear  growth 
rates  obtained  from  the  dispersion  relation,  in  Eq.  (30),  (solid  curve)  over  the  frequency  spec¬ 
trum.  These  two  independent  calculations  of  the  linear  growth  rate  are  in  excellent  agreement. 
Figure  (4)  also  compares  the  efficiency  at  saturation  obtained  by  solving  the  FEL  equations 
(circle.s  (O))  with  the  calculated  values  of  efficiency  using  electron  trapping  arguments  (dotted 
curve)  given  in  section  V.  Using  the  value  of  efficiency  for  maximum  linear  growth  rate,  we 
find  from  Eq.  (41)  that  the  saturated  vector  potential  amplitude  is  »  28  volts,  correspond¬ 
ing  to  an  efficiency  of  0.37%,  whereas  Fig.  (3)  gives  a  value  of  33  volts  for  A^,  corresponding 
to  an  efficiency  of  0.52%.  The  higher  calculated  efficiency  can  be  explained  by  the  slight  in¬ 
crease  of  the  wave  number  of  the  scattered  radiation,  just  before  saturation,  solid  curve  in 
Fig.  (S).  When  k+  increases,  the  phase  velocity  of  the  ponderomotive  potential. 
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4  (Itcreaaes.  As  the  electrons  become  trapped  at  the  bottom  of  the  potential 

well,  the  ponderomotive  wave  slows  down  slightly,  hence,  the  particles  are  able  to  transfer 
more  kinetic  energy  to  the  scattered  radiation.  This  is  clearly  a  non-linear  effect,  which  linear 
theory  could  not  predict.  The  dotted  curve  in  Fig.  (S)  is  the  variation  of  the  wave  number  of 
the  space  charge  wave.  The  effects  of  the  space  charge  wave  is  negligible,  since  in  the  strong 
pump  limit,  example  I,  the  ponderomotive  potential  is  much  larger  than  the  space  charge  po¬ 
tential  as  can  be  seen  in  Fig.  (6). 

To  understand  the  phenomenon  of  trapping,  phase  space  plots  are  a  revealing  tool.  Fig¬ 
ures  (7a-d),  are  plots  of  the  relative  time  the  particles  in  one  beamlet  cross  the  following  axial 
positions:  2  <-  0,0  m,  2.0  m,  4  0  tn,  4.3  m  and  4.S  m.  Twenty  particles  are  labeled  within  the 
beamlet.  At  the  initial  position,  z  -■  0,  the  particles  enter  at  equal  intervals  in  time  since  they 
have  uniform  axial  velocity,  v„.  At  z  -  2  m  down  stream  into  the  interaction  region,  the  parti¬ 
cles  are  in  the  linear  regime  where  the  growth  rate  of  the  scattered  radiation  is  constant.  Some 
particles  have  gained  energy  while  others  have  lost  energy  depending  on  their  phase  relation 
with  the  ponderomotive  potential.  At  z  ■■  4  m.  Fig.  (7b),  the  phase  space  plot  begins  to  show 
the  signs  of  trapping.  Many  of  the  particles  are  crossing  the  z  »  4  m  plane  at  about  the  same 
time.  However,  their  velocity  spread  is  large.  Figure  (7c),  at  z  -  4.3  m,  depicts  the  particles 
before  saturation  and  shows  definite  signs  of  trapping.  At  z  »  4.S  m,  particles  labelled  4-9  in 
Pig.  (7d)  shows  spatial  bunching  and  small  velocity  spread;  these  particles  are  deeply  trapped. 
If  the  amplitude  and  period  of  the  magnetic  pump  held  is  held  fixed,  the  scattered  radiation  will 
reach  its  maximum  value  at  this  axial  position. 

Our  non-linear  formulation  is  also  applied  to  a  case  where  the  outppt  radiation  is  in  the 
submillimeter  regime,  example  2  in  Table  1.  The  pump  wavelength  and  pump  magnetic  field 
amplitude  are  2  cm  and  2.S  kG,  respectively.  The  electron  beam  energy  is  2.6  MeV,  (yg  -  6); 
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the  beam  current  is  S  kA  and  the  beam  radius  is  0.3  cm.  The  transveme  equilibrium  velocity  is 
Vgi  <-  0.078  c,  and  the  critical  transverse  velocity  is  0.22  c.  In  the  example  we  are  bare¬ 
ly  in  the  weak-pump  retime,  since,  is  less  than  three  times  that  or  Space  charge 
effects  are,  therefore,  important  in  this  example.  Figure  (8)  is  a  plot  of  the  space  charge  and 
ponderomotive  potential  for  w  S.0.5  x  10‘^  sec'*  (X  338  fim).  This  figure  shows  that  col¬ 
lective  effects  are  of  the  same  order  of  magnitude  as  the  ponderomotive  forces. 

Figure  (9)  shows  the  amplitude  of  the  vector  potential  amplitude  of  the  scattered  radia¬ 
tion,  4(2),  and  the  spatial  growth  rate,  T  -  d(ln  4  (2))/d2,  as  a  function  of  z  for 
w  -  5.05  X  10**  sec-'. 

Comparing  the  linear  spatial  growth  rate  obtained  from  the  dispersion  equation  in  (30), 
(solid  curve  in  Fig.  (10)),  with  the  growth  rate  from  the  linear  regime  of  the  non-linear  calcula¬ 
tion  (crass  (x)),  we  again  obtain  excellent  agreement.  The  theoretical  efficiency  based  on  Eqs. 
(38)  and  (31)  (dotted  curve  in  Fig.  (10))  as  compared  with  the  results  using  the  non-linear  for¬ 
mulation,  (circles  (O)),  is  remarkably  good.  The  changes  in  wave  number  of  the  scattered  radi¬ 
ation,  ki.{z)  near  saturation  (solid  curve  in  Fig.  (11))  did  not  enhance  the  efficiency  because 
the  effect  is  balanced  by  the  increase  in  the  space  charge  potential  wave. 

The  particle  phase  space  plots  of  Figs.  (12a-c),  are  very  similar  in  nature  to  those  in  Figs. 
(7a-d).  Figure  (12a)  contains  phase  plots  at  z  -■  0.0  m  and  0.35  m  corresponding  to  the  initial 
position  and  a  point  in  the  linear  interaction  regime.  At  the  z  0.7  m  plane  just  before  satura¬ 
tion,  Pig.  (12b),  shows  the  beginning  of  particle  trapping.  Figure  (12c)  contains  the  phase 
space  plot  when  the  radiation  field  has  saturated,  z  0.77  m. 

Figure  (13)  shows  the  scaling  of  the  linear  spatial  growth  rate  and  maximum  efficiency  as 
a  function  of  the  pump  magnetic  field  amplitude.  Do,  at  a  fixed  output  frequency.  The  output 


SPRANOLE,  TANQ,  AND  MANHEtMER 


radiation  frequency  is  held  constant  by  requirini  that  y^o  Rttd  iiro  be  kept  fixed.  The  electron 
beam  and  magnetic  pump  parameters  are  basically  the  same  u  those  of  example  1  in  Table  1, 
except  that  the  magnetic  pump  amplitude  ranges  from  0.2S  kO  to  6  kO.  To  keep  the  frequency 
fixed,  while  Hq  is  varied,  the  electron  beam  energy  is  changed  such  that  is  held  at  the  con* 
stant  value  of  100.  The  output  frequency  used  for  Fig.  (13)  is  chosen  at  the  maximum  growth 
rate,  which  is  very  close  to  w  -  2y,^oc  Icq  ■>  2.S2S  x  10'^  sec~*  corresponding  to  a  wavelength  of 
k  *  0.75  ftm.  The  critical  transverse  velocity,  as  discussed  in  section  V,  occurs  for  these 
parameters  at  e  pump  magnetic  field  of  Bo  •  I.IS  kO.  Above  this  value  of  pump  field  the  PEL 
process  is  in  the  strong  pump  regime,  while  sufficiently  below  Bq  »  I.IS  kO  the  scattering  pro¬ 
cess  is  in  the  weak  pump  regime.  In  Fig.  (13)  the  crosses  (x)  denote  the  linear  spatial  growth 
rate  obtained  from  the  non-linear  simulations,  while  the  solid  curve  is  obtained  from  the 
dispersion  relation  in  Eq.  (30).  Also  in  this  figure  is  a  comparison  of  efficiency  estimated  from 
Eq.  (38)  using  trapping  arguments  (dashed  curve)  and  actual  numerical  simulation  results  (cir¬ 
cles  (O)). 

b)  Efficiency  Enhancement  by  Contouring  Magnetic  Pump  Period 

According  to  Eqs.  (24)  and  (26)  the  phase  velocity  of  the  total  longitudinal  wave  poten¬ 
tial,  i.e.  ponderomotive  plus  space  charge  is 

-  <u/(k+  +  *,)  (42) 

where  u  and  Ac 4.  are  the  radiation  frequency  and  wave  number  and  k„  -  2v/l  is  the  wave 

number  of  the  pump  field.  It  has  been  assumed  in  writing  (42)  that  the  wavenumber  of  the 
ponderomotive  and  space  charge  waves  are  identical.  The  longitudinal  wave  potential  is  respon¬ 
sible  for  axially  bunching  and  eventually  trapping  the  electrons.  If  the  magnetic  pump  period  is 
held  fixed,  the  radiation  field  reaches  its  maximum  value  when  the  electrons  are  trapped  at  the 
bottom  of  the  longitudinal  potential  wells,  as  can  be  seen  for  example  in  Fig.  (7d).  Just  before 
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the  radiation  field  saturates,  the  electroiu:  are  somewhat  spatially  bunched  and  trapped  near  the 
bottom  cf  the  wave  potential,  see  Pig.  (7c).  The  trapped  electrons  at  this  point  can  be  con¬ 
sidered,  for  our  purpose,  to  form  a  macro-particle.  By  appropriately  reducing  the  phase  velocity 
in  Eq.  (42)  as  a  function  of  axial  distance  down  the  interaction  region,  the  kinetic  energy  of 
this  niacro-particle  can  be  further  reduced  and  converted  into  wave  energy.  The  phase  velocity- 
must  be  reduced  in  such  a  way  so  that  the  inertial  potential  of  the  trapped  macro-particle  is  al¬ 
ways  less  than  the  potential  of  the  growing  longitudinal  wave.  According  to  Eq.  (42),  the  phase 
velocity  can  be  reduced  by  decreasing  the  period  of  the  magnetic  pump  as  a  function  of  z.  In 
order  for  the  macro-particle  to  remain  trapped,  the  spatial  r<.te  of  change  of  the  pump  period 
must  be  sufficiently  slow.  In  principle  virtually  all  the  kinetic  energy  of  the  macro  ,jarticle  can 
be  extracted  and  converted  to  wave  energy.  However,  not  all  the  beam  particles  comprise  the 
macro-particle;  some  are  untrapped.  Converting  particle  kinetic  energy  into  radiation  by  vary¬ 
ing  the  wave  velocity  is  somewhat  analogous  to  the  reverse  process  of  particle  acceleration  in 
say  an  RF  linac.  In  a  wave  accelerator,  the  energy  associated  with  the  accelerating  slow  elec¬ 
tromagnetic  wave  is  converted  into  particle  kinetic  energy.  However,  the  wave  energy  in  these 
accelerators  does  not  decay,  since  it  is  continuously  resupplied  by  external  microwave  sources. 

We  will  illustrate  efficiency  enhancement  by  contouring  the  pump  period  while  holding 
the  amplitude  of  the  pump  magnetic  vector  potential  constant,  using  the  parameters  of  example 
1  in  Table  1.  The  same  principle  of  efficiency  enhancement  can  also  be  applied  to  example  2. 
Figure  (7c)  shows  that  a  z  4.3  m,  the  electrons  are  somewhat  spatially  bunched  at  the  optical 

wavelength  X  ~  0.75Mm  and  the  radiation  field  is  nearly  saturated.  At  this  point,  we  simply  in¬ 

creased  the  pump  wavenumber  koiz)  exponentially  as  a  function  of  z  instead  of  optimally  con¬ 
touring  the  pump  period  in  units  of  cm~'  according  to  the  empirical  formula. 

2.r 

1.5  «  <  z* 

+  [exp(0.002(r-z*))-ll  ^  ^ 
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where  Zk  ^  4.3  m.  The  period  of  the  magnetic  pump,  /(z),  is  depicted  in  Fig.  (14).  The  spatial 
decrease  of  /  results  in  a  large  increase  in  the  amplitude  of  the  wave  vector  potential  as  shown 
in  Fig.  (14).  For  this  particular  case,  the  contouring  is  terminated  at  z  ■■  13  m  and  the 
efficiency  at  this  point  is  already  20%.  In  principle,  the  pump  wavelength  contouring  can  be 
continued  and  even  higher  efficiencies  achieved.  Figures  (ISa-d)  are  the  phase  plots  with  con¬ 
touring  at  2  —  5  m,  7  m,  10  m  and  13  m.  At  z  -  5  m,  the  majority  of  the  particles  are  well 
bunched.  At  z  —  7  m,  12  out  of  20  particles  are  trapped  by  the  ponderomotive  potential  wells; 
the  same  12  particles,  remain  trapped  even  at  z  -■  13  m.  Since  the  amplitude  of  the  pondero¬ 
motive  potential  is  proportional  to  the  radiation  held  it  increases  as  the  radiation  held  increases. 
Once  the  particles  are  trapped  the  particles  remain  trapped  and  continually  lose  energy  if  the 
pump  period,  /,  is  decreased  adiabatically. 

It  should  be  noted  that  it  is  not  appropriate  to  simply  increase  the  pump  magnetic  held 
amplitude  as  a  function  of  axial  position  in  order  to  enhance  efficiency.^''  Increasing  spatially 
deepens  the  ponderomotive  potential  well  but  it  also  slows  down  the  axial  electron  velocity  and 
thus  synchronism  is  lost.  To  maintain  synchronism  ko  must  also  be  increased. 
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Appendix  A 

GENERAL  NON-LINEAR  FORMULATION 

In  this  appendix,  we  outline  the  general  formulation  of  the  PEL  equations  taking  into  ac¬ 
count  spatial  harmonics  in  the  magnetic  pump  field  as  well  as  spatial  and  temporal  harmonics  in 
the  scattered  fields.  Furthermore,  the  polarization  of  the  electromagnetic  field  is  arbitrary  and 
permitted  to  evolve  according  to  the  non-linear  particle  field  dynamics. 

The  vector  potential  of  the  periodic  pump  field  containing  spatial  harmonics  of  variable 
amplitudes  and  wave  numbers  is  expressed  as 

A<,(z)  -  I  cosjiM  +  sin  |/m  J  (Al) 

where  the  amplitude  and  fundamental  wavenumber  are  slowly  varying  function  of  2.  This  field 
is  not  curl  free,  but  is  a  gbod  approximation  to  the  exact  helically  symmetric  field  near  the 
r  —  0  axis,  when  rnkgr^  <  1,  where  is  the  radius  of  the  electron  beam.  The  pump  magnetic 
field  associated  with  (Al)  is  given  by 


where 


and 


*0^2)  -  X  I  cosjm  kjz')dz'  +  (p„(2) j 

+  sio|m  J^ko(z’)dz'  +  ^;„(z)j  ^^1, 


B,.Jz) 


(nt  ko{z)Ao^„{z))'^  + 


I  dz 


ifi„(z)  -tan  ' 


9A„„{z) 


Bz 


(m  k„(z)Ao„(z)) 
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are  slowly  varying  functions  of  z.  The  period  of  the  mth  spatial  harmonic  of  the  pump  magnetic 
field  is  a  function  of  z  and  is 


where 


/m(2) 


m  k„iz)  + 


dz 


m  ke(z)  » 


dz 


Similarly  the  general  form  for  the  scattered  electromagnetic  field  and  electrostatic  field  in 
terms  of  the  vector  potential  A  (z,r),  and  scalar  potential  (z.t) 


A  (z,t)  -  ^  ~ 

/l"l  \ 

+  Ay,„(z)  sin  I  n  f^ky„(.z')dz'  —  ntat  +  e^j. 


and 


(A2) 


where  the  amplitudes  of  the  potentials  Ax„(.z),  Ay„iz),  and  ^t(z)  as  well  as  wave  numbers 
kx,„{z),  ky  „(z)  and  kiiz)  are  slowly  varying  functions  of  r. 


Using  the  same  procedure  as  used  to  derive  Eqs.  (7a-f)  we  obtain  the  following  set  of 
equations  for  the  spatially  slowly  varying  amplitudes  and  wavenumbers. 

~  Ax,„U)  -  ~“X”  cos|fiX*^,rt'^2')*'  “  +  ^x.n\dto- 


In  kl'^iz)  i4,,„(z)Ac;,‘iHz)|  -  sinjnX  ky,„{z')dz'  -  nwr  +  9^„^dto 


Ay  „{z)  -  JyizA)  sin|rtX,*j',/i(2')<fe'  -  n<ot  + 
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-  -—X*  cos(«XV,(z')*'  -  nmt  +  9y  „]dt,. 

-  4^^“  </,(z.<)  sin J/ X  ~ 

/*  ki  ^,(z)  -  4  “  Ji(2.r)  cosj/X,  k,iz’)d2'  -  /w»  +  (A4a-0 

Subsiituting  the  expressions  for  i(z,i)  from  Eq.  (13),  into  the  Eqs.  (A4a-0  and  integrating,  we 


obtain 


-  «C,^„(2)  .4x,^(2) 


- ^Iniz)  Ay,„(z) 


.  -4UI-,.  -  n,.rU,.z)  +  »,,.lrfv 

2«  -^1 

«  _4l^|„  cosfrt  +  »y,n\dto 

’  ’  c  do  HzUo>rUo,z))  * 


/  .  -4|elno  Vjo  X  sin|/X^  Mz')*'  ~  /‘«t(/o,z)  +  0z,/j<f/o- 

/^/c/  <A,(r)  -  -4|c|«„  v,o  sr  cos  [//\(z')«/2'  -  /a)T(r„,z)  +  (ASa-O 


#  •  ’  :*«>.  ’Hi'  .''"  v 

.  v.; 
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The  expression  of  the  particle  orbit  equation,  (14),  remains  the  same.  Equations  (14), 
(ASa-f).,  and  along  with  the  definitions  of  momenta, 

7j,(2,t)  -  (z)  +  Axiz.r)), 

t)j,(z.t)  —  -^(Xo.)-  (2)  +  Ayiz.r)), 

and 

T)i(2,  t)  -  moyiz,  t)  Kj (2,  t)  ,  (A6) 

form  the  full  set  of  self-consistent  PEL  equations. 

Setting  M  -  m  —  /  -  1,  and  requiring  *,,1(2)  -  *^.1(2)  —  k+iz)  and  Ax,\(z)  —  Ay_\(z)  — 
Aiz),  equations  (A5a-f)  reduce  to  the  fundamental  harmonic  equations  in  (22a, b)  and  (23a,b) 
for  a  circularly  polarized  electromagnetic  field. 
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Table  I  —  Optical  and  Submillimeter  Illustrations  of  FELs 
(Constant  Magnetic  Pump  Parameters) 


Magnetic  Pump  Parameters 

Example  #1 

Example  #2 

Pump  Wavelength 

/ 

l.S  cm 

2.0  cm 

Pump  Amplitude 

Bo 

6,0  kG 

2.5  kG 

Electron  Beam  Parameters 

Beam  Energy 

Eo 

66  MeV  (y,  -  131) 

2.6  MeV(y„  -  6) 

Beam  Current 

lo 

2  kA 

5  kA 

Axial  Gamma 

y«» 

100 

5.4 

Beam  Radius 

^0 

0.1  cm 

0.3  cm 

Equil.  1  Velocity 

BoL 

6.4  X  10"^ 

0.078 

Critical  1  Velocity 

Bern 

l.S  X  10-^ 

0.22 

Beam  Strength  Parameter 

0.14 

0.87 

Self  Potential  Energy  Spread 

BElEo 

0.08% 

1.7% 

Output  Radiation  Parameters 

Radiation  Wavelength 

X 

0.75  Min 

338  Min 

Linear  enfolding  length* 

% 

1 

1 

3 

38  cm 

5.3  cm 

Efficiency* 

V 

0.52% 

9.2% 

Saturated  A-Fieid* 

A 

33  volts 

7.4  X  10^  volts 

Radiation  Power* 

Po 

0.69  GW 

1.2  GW 

*For  maximum  growth  rate. 
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Pi|.  1  —  Sciwiiittic  of  the  fTeo-electton  User  model.  The  unmodulated 
electron  beam  entera  the  Interaction  region  from  the  leit.  The  pump 
AeM  builda  up  adiabatkatly  and  readiea  a  conatant  amplitude  for  t  >  0. 


Pig.  2  —  .A  diagram  to  illvwirate  the  aymmetry  property  of  the 
function  GIiq.  z,  t)  in  the  (r.  to)  plane. 
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Fit-  3  —  Wtve  vector  potential  AU)  and  spatial  linear  trowth  rate  F  as  a  ftinction  of 
axL^I  distance  for  example  1  in  the  optical  regime.  The  frequency  is  chosen  to  give 
the  maximum  linear  spatial  growth  rate. 


Fig.  4  —  A  comparison  of  the  growth  rate  in  the  linear  regime  of  the  non-linear  simulation 
(crosses  (x))  with  the  growth  rate  from  linear  theory  (solid  curve),  and  a  comparison  of 
efficiency  from  non-linear  theory  (circles  (O))  with  that  from  linear  theory  using  trapping 
arguments  (dashed  ctrve)  as  a  function  of  frequency  for  example  I. 
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1.S  3.0  4.5  6.0 

AXIAL  DISTANCE  (METERS) 


Fi(.  S  -  The  verleilon  in  wavenumber  of  the  scetteied  redietlon  k^(t)  end  the 
variation  in  wavenumber  of  the  apace  charte  potential  kU)  u  a  function  of 
axial  diatanoa  for  example  I  at  the  frequency  correapondint  to  maximum  linear 
spatial  growth. 


AXIAL  DISTANCE  (METERS) 

Fig.  6  —  A  comparison  of  the  magnitude  of  the  ponderorootive  potential  |4^(r)|  and 
the  space  charge  potential  |4(r)|  as  a  function  of  axial  distance  for  example  I  at  the  fre¬ 
quency  corresponding  to  msximum  iinoar  spatial  growth. 


35 


MZ)  (10  ®  VOLTS) 


(oas/uio  .,01)  -  (2)*  A 


NRL  MEMORANDUM  REPORT  4034 


(b) 

'ig.  7  (Continued)  —  Phase  space  plots  of  velocity  versus  the  relative  time  the  particles  in  one  beamlet  cross 
the  following  axial  positions  (a)  r  ••  0.0  m  and  z  —  2.0  m,  (b)  z  ••  4.0  m,  (c)  z  -  4.3  m  and  (d)  z  45  m  for 
example  < 
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WtUo.z)  -  u>t(o,z) 

(c) 

Fig.  7  (Continued)  —  Phase  space  plots  of  velocity  versus  the  relative  time  the  particles  in  one  beamlet  cross 
the  following  axial  positions  (a)  z  —  0.0  m  and  z  -  2.0  m,  (b)  z  4.0  m,  (c)  z  ••  4.3  m  and  (d)  z  —  4S  m  for 
example  1 
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a»T(to,z)  a)T{o,z) 

(d) 

Fig.  7  (Continued)  —  Piiase  space  plots  of  velocity  versus  the  relative  time  the  particles  in  one  beamlet  cross 
the  following  axial  position:  (a)  z  0.0  m  and  z  —  2.0  in,  (b)  z  -  4.0  m,  (c)  z  —  4.3  m  and  (d)  z  -  45  m  for 
example  1  ^ 

*^4 
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Fig.  8  —  A  comparison  of  the  magnitude  of  the  ponderomotive  potential 

the  space  charge  potential  |4i(z)|  as  a  function  of  axial  distance  for  example  2  at  the 

quency  corresponding  to  maximum  linear  spatial  growth. 


Fig. '  —  The  wave  vector  potential  A  (z)  and  the  spatial  growth  rate  F  as  a  lunction  of 
axial  distance  for  example  2  at  the  frequency  corresponding  to  maxinum  I' near  growth. 
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Fig.  10  —  A  comparison  of  the  growth  rate  in  the  linear  regime  of  the  non-linear  simula¬ 
tion  (crosses  (x))  with  the  growth  rate  from  linear  theory  (solio  curve),  and  a  compari¬ 
son  of  efficiency  from  the  non-linear  theory  (circles  (O))  with  that  from  linear  theory 
using  trapping  arguments  (dashed  curve)  as  a  function  of  frequency  for  example  2. 


Fig.  11  —  The  variation  in  wavenumber  of  tho  scattered  radiation  k^.{z)  and 
the  variation  in  wavonnmber  of  the  space  charge  potential  f.  (z)  as  a  function  of 
axial  distance  for  example  2  at  ths  frequency  corresponding  to  maximum  linear 
spatial  growth. 
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Fig.  12  —  Phase  space  plots  of  velocity  versus  the  relative  time  the  particles  in  one  beamlet  cross  the  following 
axial  positions  (a)  r  -  0.0  m  and  r  —  0.35  m,  (b)  z  -  0.7  m  and  (c)  z  —  0.77  m  for  example  2 
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Fig.  12  (Continued)  —  Phue  spice  plots  of  velocity  versus  the  relitive  time  the  particles  in  one  beimlet  cross  the 
following  ixiil  positions  (•)  z  -  0.0  m  and  z  —  0.3S  m.  fb)  z  •  0.7  m  and  (c)  z  —  0.77  m  for  example  2 
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FI*  13  —  A  comparison  of  the  maximum  trowth  rate  in  the  linear  re|ime  of  the  non¬ 
linear  simulation  (crosses  (x))  with  the  maximum  prowth  rate  Oom  linear  theory  (solid 
curve),  and  a  comparison  of  efflcieacy  from  non-linear  theory  (circles  (O))  with  that 
from  linear  theory  using  i»e  trapping  arguments  (dsshed  curve)  as  a  (Vinction  of  the  mag¬ 
netic  pump  held  amplitude  for  a  fixed  output  Oequency. 


AXIAL  DISTANCE  (METERS) 


Fig.  14  —  Enhancement  of  radiation  field  by  decreasing  the  magnetic  pump 
period.  The  eiBciencv  has  increased  from  0.52%  at  z  -  4.5  m  with  a  constant 
pump  period  to  30%  at  i  13  m  with  the  period  of  the  pump  changing  as 
shown. 
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Fi»- 15  (Continued)  -  Pheee  plots  with  contoured  ratcnetic  pump  period  shown  in  Fip.  (14)  ut  (•) 

(c)  *  -  10  m  and  (d)  s  -  13  m 
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